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Abstract 



We introduce a new dynamical indicator of stability based on the Extreme Value statistics 
• • showing that it provides an insight on the local stability properties of dynamical systems. The 

indicator perform faster than other based on the iteration of the tangent map since it requires 
only the evolution of the original systems and, in the chaotic regions, gives further information 
about the information dimension of the attractor. A numerical validation of the method is 
presented through the analysis of the motions in a Standard map. 

1 Introduction 

The analysis of stability for discrete and continuous time dynamical systems is of fundamental im- 
portance to get insights in the dynamical structure of a system. The distinction among regular and 
chaotic orbits can be easily made in a dissipative case whereas for conservative systems it is usually 
an hard task especially when we are dealing with many degrees of freedom or when the phase space 



is splitted into chaotic regions and regular islands. A large number of tools known as indicators of 
stability have been developed to accomplish this task: Lyapunov Characteristic Exponents (LCEs) 



(Wolf et al., 1985), (Rosenstein et al. 1993), (Skokos 2010) and the indicators related to the Return 
Time Statistics flKac[ |1934| ), ( |Gao| |1999| ), flHu et al.\ |2004| ) , ( |Buric et al] |2005[ ) have been used 
from a long time as efficient indicators. Nevertheless, in the recent past, the need for computing 
stability properties with faster algorithms and for systems with many degrees of freedom resulted in 
a renewed interest in the technique and different dynamical indicators have been introduced. The 



Smaller Alignment Index (SALI) described in Skokos et al. (2002) and Skokos et al. (2004) ,the 



Generalized Alignment Index (GALI), introduced in Skokos et al. (2007 land the Mean Exponential 



Growth factor of Nearby Orbits (MEGNO) discussed in Cincotta et al. ( 2003 ) Gozdziewski et al. 
(2001) are suitable to analyse the properties of a single orbit. They are based on the divergence 



of nearby trajectories and require in principle the knowledge of the exact dynamics. Another class 
of indicators is based on the round off error properties and has been discussed in |Faranda et al. 
(2011c): the divergence between two trajectories starting from the same initial condition but com- 
puted with different numeric precision can be used to overlight the dynamical structure. The so 
called Reversibility Error that measure the distance between a certain initial condition and the end 
point of a trajectory iterated forward and backward for the same number of time steps give basically 
the same informations. 
Other tools such as Fidelity and Correlation decay can be successfully used to characterize stability 



properties of ensemble of orbits as explained in Liverani et al. (2007) and Turchetti et al. (2010) 



Eventually, Frequency Map Analysis has been used to provide informations about the resonance 



structure of a system (Laskar 1999), (Robutel & Laskar, 2001). 



Although these indicators have been developed to accomplish the same task, each of them 
presents specific features and limitations and often it is necessary to combine a certain number of 
indicators to get quantitative information about the dynamics. Another aspect to consider when 
comparing the efficiency of chaos indicators is the computational cost: each variational method 
(SALI, MEGNO, mLCE) needs to iterate both the map and the tangent map forward during n 
steps. The latter map, although linear, is the computationally most expensive due to the fact that 
it needs the evaluation of the Jacobian matrix at every step. On the other hand the Round off error 
methods are the computationally less expensive as they require only the iteration of the dynamics. 
Fidelity and Correlations are usually computed using Monte Carlo simulations and therefore they 
are inaccessible for systems with many degrees of freedom Turchetti et al. (2010). Therefore, it is 



clear that there is still the need to introduce versatile indicators that distinguish regular from chaotic 
behaviors and possibly give further informations on the dynamics. The pourpose of this paper is 
to use the series of extrema of some specific observable computed using the orbits of dynamical 
systems as a new indicator of stability. 



Extreme Value Theory was originally introduced by Fisher &; Tippett (1928), Gnedenko (1943)to 



study the maxima of a series of independent and identical distributed variables: under very general 
hypothesis a limiting distribution called Generalised Extreme Value (GEV) distribution exists for 
the series of extremes. An extensive account of recent results and relevant applications is given in 



Ghil et al. (2011 ). In the recent past this theory has been adapted to study the output of dynamical 



systems. As we will explain in detail in the next section it is not trivial to observe the asymptotic 



GEV distribution in dynamical systems: some sort of independence of maxima must be recovered by 
requiring certain mixing conditions on the orbits. Furthermore, we need to introduce some peculiar 
observables that satisfy the condition proposed by Gnedenko on the parent distribution of data: 
they are related to the closest return of a trajectory in a ball centered around the starting point 
and therefore allow a very detailed description of the dynamics in the neighborhood of the initial 
condition. 

If all these requirement are satisfied it is possible to observe an Extreme Value statistics that 
converges to the GEV distribution family. The parameters of the distribution are intimately related 



with some relevant dynamical quantities such as the dimension of the attractor (Freitas et al. 2009) 



(Faranda et al. 2011a). We will use these features to show the reliability of GEV distribution 



parameters in discriminate regular from chaotic behaviors pointing out the further informations 
that is possible to extract regarding the dynamics. 

The paper is organised as follows: in section 2 we explain how to introduce an EV statistics in 
dynamical systems pointing out the difference between regular and chaotic orbits. In section 3 
we describe the numerical algorithm and procedure used to compute the parameters of the GEV 
distribution. Eventually, in section 4 we perform some tests on the Standard Map to validate 
numerically the use of GEV parameters as dynamical indicators. 

2 Extreme Statistics as dynamical indicator 
2.1 Extreme Value Theory in dynamical systems 



Gnedenko (1943) studied the convergence of maxima of i.i.d. variables 

Xq,Xi, ..X m -i 

with cumulative distribution (cdf ) F{x) of the form: 

F{x) = P{a m (M m - b m ) < x} 

Where a m and b m are normalizing sequences and M m = max{Vo, Xi, ..., V m _i}. It may be rewritten 
as F{u m ) = P{M m < u m } where u m = x/a m + b m . Under general hypothesis on the nature of the 



parent distribution of data, Gnedenko (1943) show that the distribution of maxima, up to an affine 



change of variable, obeys to one of the following three laws: 
• Type 1 {Gumbel). 



Type 2 (Frechet). 



Type 3 (Weibull). 



E(x) 



E(x) 



E{x) = exp(— e x ), — oo < x < oo 

Jo, x< 

I exp(— x~€), for some £ > 0, x > 

{exp(— (— x)~£), for some £ > 0, x < 
1, x>Q 



(1) 



(2) 



(3) 



Let us define the right endpoint xp of a distribution function F(x) as: 

xp = sup{x : F(x) < 1} 



(4) 



then, it is possible to compute normalizing sequences a m and b m using the following corollary 
of Gnedenko's theorem : 

Corollary (Gnedenko): The normalizing sequences a m and b m in the convergence of normalized 
maxima P{a m (M m — b m ) < x} — > F(x) may be taken (in order of increasing complexity) as: 



where 



• Type 1: a m = [G(7 m )]~ 1 , b m = -y m ; 

• Type 2: a m = 7" 1 , b m = 0; 

• Type 3: a m = (xp - 7m)" 1 , b m = x F ; 

7 m = F _1 (l - 1/rn) = inf{x; F(x) > 1 - 1/m} 



G(t) 



XF 1 - F(u) 



du, t < xp 



(5) 



(6) 



In (Faranda et al., 2011b) we have shown that this approach is equivalent to fit unnormalized 



data directly to a single family of generalized distribution called GEV distribution with cdf: 



Fg{x]^,(t,0 =exp 



i + e' 



x — IX 

a 



-W 



(7) 



which holds for l+£'(x— fi)/cr > 0, using /i£l (location parameter) and a > (scale parameter) 



as scaling constants in place of b m , and a m (Pickands III, 1968), in particular, in Faranda et al. 



(2011b) we have shown that the following relations hold: 



(i = b r , 



a 



£' E R is the shape parameter also called the tail index: when £' — > 0, the distribution cor- 
responds to a Gumbel type ( Type 1 distribution). When the index is positive, it corresponds to 
a Frechet (Type 2 distribution); when the index is negative, it corresponds to a Weibull (Type 3 
distribution). 

In the last decade many works focused on the possibility of treating time series of observables 



of deterministic dynamical system using EVT. For example, Balakrishnan et al. (1995) and more 



recently Nicolis et al. (2006) and Haiman (2003) have shown that for regular orbits of dynamical 



systems we don't expect to find convergence to EV distribution. 

The first rigorous mathematical approach to extreme value theory in dynamical systems goes back 



to the pioneer paper by P. Collet in 2001 (Collet, 2001). Important contributions have successively 



been given by Freitas & Freitas (2008), Freitas et al. (2009).Freitas et al. (2010) and by Gupta et al. 



(2009). The starting point of all these investigations was to associate to the stationary stochastic 



process given by the dynamical system, a new stationary independent sequence which enjoyed one of 
the classical three extreme value laws, and this laws could be pulled back to the original dynamical 
sequence. 

Let us consider a dynamical systems (fi, B, u, /), where Q is the invariant set in some manifold, 
usually M d , B is the Borel a-algebra, / : $7 — > 0, is a measurable map and v a probability /-invariant 
Borel measure. 

In order to adapt the extreme value theory to dynamical systems, we will consider the stationary 
stochastic process Xo,Xi, ... given by: 

X m (x) = g(dist(f m (x),()) VmeN (8) 

where 'dist' is a distance on the ambient space f2, £ is a given point and g is an observable 
function, and whose partial maximum is defined as: 

M m = max{X ,...,X m _ 1 } (9) 

The probability measure will be here the invariant measure v for the dynamical system. We will 
also suppose that our systems which verify the condition D2 and D' which will allow us to use the 
EVT for i.i.d. sequences. As we said above, we will use three types of observables gi,i = 1,2,3, 
suitable to obtain one of the three types of EV distribution for normalized maxima: 

5 i(x) = -log(dist(x,C)) (10) 

52 (x) = dist(x,C)- 1/a (11) 

g 3 (x) = C -dist(x,() 1/a (12) 

where C is a constant and a > £ M. 

Using these observables we can obtain convergence to the Type 1,2,3 distribution if one can 
prove two sufficient conditions called D2 and D' and which we briefly explain here: these conditions 
basically require a sort of independence of the stochastic dynamical sequence in terms of uniform 
mixing condition on the distribution functions. In particular condition D2, introduced in its actual 



form by Freitas-Freitas Freitas &i Freitas (2008), could be checked directly by estimating the rate 

of decay of correlations for Holder observables. 

If X m ,m > is our stochastic process, we can define Mji = {Xj,Xj + i,- ■ ■ ,Xj+i} and we put 

M Qtm = M m . 

The condition D2(u m ) holds for the sequence X m if for any integer I, t, m we have 

\v(X > u m ,M t j < Um) - v(Xq > u m )v(M t j < u m )\ < y(m,t) 

where 7(771, t) is non-increasing in t for each 777 and 7777(777, t m ) — > as 777 — > 00 for some sequence 

t m = o(m), t m -> OO. 



We say condition D'(u m ) holds for the sequence X m if 



[m/l] 



lim limsupm N v{Xq > u m ,Xj > u m ) = 

Z->oo m t - f 



Instead of checking the previous conditions, we can use another results that established a connection 
between the extreme value laws and the statistics of first return and hitting times, see the papers by 



Freitas et al. (2009) and Freitas et al. (2011 ). They showed in particular that for dynamical systems 



preserving an absolutely continuous invariant measure or a singular continuous invariant measure u, 
the existence of an exponential hitting time statistics on balls around v almost any point £ implies 
the existence of extreme value laws for one of the observables of type gi, i = 1, 2, 3 described above. 
The converse is also true, namely if we have an extreme value law which applies to the observables 
of type gi,i = 1,2,3 achieving a maximum at £, then we have exponential hitting time statistics 
to balls with center £. Recently these results have been generalized to local returns around balls 



centered at periodic points (Freitas et al. 2010). 



2.2 Extreme Value Statistics in mixing and regular maps 



In Faranda et al. (2011b) and Faranda et al. (2011a) we have analised both from an analytical and 



numerical point of view the Extreme Value distribution in a wide class of low dimensional maps 
showing that, when the conditions D' and D2 are verified, the block maxima approach can be used 
to study extrema. It consists of dividing the data series of length k of some observable into n bins 
each containing the same number m of observations, and selecting the maximum (or the minimum) 
value in each of them (Coles et al. 1999). Using gi observable functions we have shown that a 



first order approximation of the GEV parameters in mixing maps can be written in terms of m (or 
equivalently n) and the dimension of the attractor d: 
For gi type observable: 



a = - a hxik/n) f ' = 

a a 



For 52 type observable: 



For 53 type observable: 



a ~ n' 1 '^ 11 ~ n- 1 /^ £' 



ad 



a ~ n 



l/(ad) 



H = C £' 



ad 



(13) 



(14) 



(15) 



while the higher order terms contain explicit dependence on the density measure. 



For regular maps independently on the observable chosen, for periodic or quasi-periodic orbit 



we do not observe convergence to the GEV distribution. In Faranda et al. (2011b) we have analysed 



what happens when observable gi are used to study empirical maxima distribution in regular maps: 



Nicolis et al. (2006) have shown how it is possible to obtain an analytical EV distribution which 



does not belong to GEV family choosing a simple observable: they considered the series of distances 
between the iterated trajectory and the initial condition: 

Y m (x = ft) = dist(/ 4 C, M m = min{Y , ...Y m _i} 
For this observable they have shown that the cumulative distribution F(x) = P{a m (M m — b m ) < 



x} of a uniform quasi periodic motions is not smooth but piecewise linear (Nicolis et al. (2006), 
Figure 3). F(x) must correspond to a density distribution continuous obtained as a composition of 
Heaviside step functions: each box must be related to a change in the slope of F(x). 
By applying the observable gi we just remodulate this piecewise linear F{x) but it is clear that 
we don't obtain any kind of convergence to the GEV distribution. In terms of density functions, 
we can observe multimodal distributions and the number of modes and their positions are highly 



dependent on both n,m and initial conditions Faranda et al. (2011b) 



Using Nicolis et al. (2006) results it is possible to understand the empirical distribution obtained 
once applied gi functions: since density distribution of M m is a composition of box functions, when 
we apply gi observables we modulate it changing the shape of the boxes. Therefore, we obtain a 
multi modal distribution modified according to the observable functions gi. In the case of pure 
periodic motion this kind of extreme value distribution must asymptotically be a Dirac delta as we 
pick up always the same Yj. 



3 The numerical algorithm 

As we have already said to observe a GEV distribution of maxima orbits must satisfie D' and 
Z?2 conditions or an exponential decay of the Hitting Time Statistics whereas for periodic or quasi- 
periodic motions we have different distributions but never a GEV. This give us a way to discriminate 
the kind of motion simply looking at the extreme value statistics. From a practical point of view 
we can introduce a simple algorithm to perform this task: 

1. Compute the orbit of the dynamical system for k iterations. 

2. Compute the series X m (x) = g(d\st(f m {x), ()) where £ is the initial condition. 

3. Divide the series in n bins each containing m data. 

4. Take the maximum in each bin and test if this empirical distribution give the parameters 
expected by the theory. 

In the next section we will describe how to use it operationally in a meaningful example of a 
dynamical systems that present coexistence of regular and chaotic motions and therefore allow to 
test the validity of our indicator: the Standard map. Before presenting the results, we need to 
clarify the numerical inference procedure that we use to obtain the parameters of the GEV distribu- 



tion. In Faranda et al. (2011b) we have used a Maximum Likelihood Estimation (MLE) procedure 
working both on pdf and cdf (cumulative distribution function), since our distributions were ab- 
solutely continuous and the minimization procedure was well defined. In a general case, when we 
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are dealing also with singular measures, we may have a cdf which is not anymore absolutely con- 



tinuous and consequently the fitting procedure via MLE could give wrong results (Faranda et al. 



2011a). To avoid these problems we have used an L-moments estimation as detailed in Hosking 



(1990). This procedure is completely discrete and can be used both for absolutely continuous or 



singular continuous cdf. The L-moments are summary statistics for probability distributions and 
data samples. They are analogous to ordinary moments which meant that they provide measures 
of location, dispersion, skewness, kurtosis, but are computed from linear combinations of the data 
values, arranged in increasing order (hence the prefix L). Asymptotic approximations to sampling 



distributions are better for L-moments than for ordinary moments Hosking (1990), Figure 4]. The 



relationship between the moments and the parameters of the GEV distribution are described in 



Hosking (1990), while the 95% confidence intervals has been derived using a bootstrap procedure. 



We have explained in Faranda et al. (2011a) the issues of applying directly a Kolmogorov Smirnov 



or a Chi-square test to see if the data really belongs to the GEV distribution when the measure is 
singular: even if both tests fail the GEV model is reasonable as it is the closest continuous repre- 
sentation of the empirical discrete distribution we get. 

Therefore, instead of using as a dynamical indicator the goodness of the fit to the continuous GEV 
model, we will check the deviations of the parameters with respect to the theoretical expected values. 

Our results will be studied against the Divergence of two nearby trajectories and the Reversibil- 
ity error. In Faranda et al. (2011c) we have shown that these indicators give insights into the 



structure of a dynamical system as well as others such as SALI, MEGNO and mLCE with which 
these indicators have been compared. 



We briefly recall here the definitions and we remind to Faranda et al. (2011c) for further clarifica- 
tions. 

The arithmetic operations such as sums or multiplications imply a round off, which propagates 
the error affecting each number. Round off algebraic procedures are hardware dependent as detailed 
in Knuth (1973). Unlike the case of stochastic perturbations, the error strongly depends on x. 
Suppose we are given a map f t (x) then we will indicate with fl(x) the correspondent numerical 



map both evaluated at the t-th iteration. The Divergence of orbits is defined as: 

A t = dist(/£(z),/i>(x)), 



(16) 



where f% and f l D stand for single and double precision iterations respectively, and dist is a suitable 
metrics. 

If the map is invertible we can also define the Reversibility error as 



R t = dist(f- t ofl( x ),x) (17) 

which is non zero since the numerical inverse f~ l of the map is not exactly the inverse of /* 
namely f~ l o f*{x) ^ x. Obviously the reversibility error is much easier to compute than the diver- 
gence of orbits (if we know explicitly the inverse map) and the information it provides is basically 
the same as the latter. Both quantities give an average linear growth for a regular map together 
with an exponential growth for a chaotic map having positive Lyapounov exponents and strong 
mixing properties. When computing Rt we will set /* = fs in order to compare with At. 



4 A case of study: the Standard Map 

The Standard map (also known as Chirikov- Taylor map or Chirikov standard map) is an area- 
preserving chaotic map defined on the bidimensional torus. It can be thought as a stick that is free 
of the gravitational force, which can rotate frictionless in a plane around an axis located in one of 
its tips, and which is periodically kicked on the other tip. This mechanical system is usually called 
a kicked rotator. It is defined by: 

f Vt+i =Vt~ |f sin(2vrx t ) mod 1 
\x t+ i = xt + yt + l mod 1 

Standard map is one of the most widely-studied examples of dynamical chaos in physics. It can 
be regular or chaotic, depending on the strength of the impulses: stronger kicks lead to chaotic be- 
haviors. The variables y and x respectively represent the angular position and angular momentum 
of the stick at the t-th kick. 

For K « 1 the motion follows quasi periodic orbits for all initial conditions, whereas if K » 1 
the motion turns to be chaotic and irregular. An interesting behavior is achieved when K ~ 1: in 
this case we have coexistence of regular and chaotic motions depending on the initial conditions 
chosen. 

First of all we have taken an ensemble of 500 initial conditions centered around Xo,yo m a small 
subset of the bidimensional torus. After iterating the map for k = 10 6 iterations, we have selected 
n = 1000 maxima each inm = 1000 observations. These values are in agreement with what emerged 



from the numerical study presented in Faranda et al. (2011b). Once computed the parameters of 



GEV distribution for each realization, we have averaged them over the different initial conditions 
checking their convergence towards the expected theoretical values when K is varied from K = 10~ 4 
up to K = 10 2 . As indicators, looking at the equations 13||l5 we have selected the shape parameters 



for the three type observables £,' (gi) , £' (g2) , £,' (93) and the scale parameter for the type 1 observable 
cr(gi). From all the parameters set, these are the ones that have a dependence on the dimension of 
the attractor d but do not have a dependence of m, therefore, once we are in the asymptotic regime, 
the results are independent on the number of observations in each bin. Results are presented in 
figure [T] In this example we set a = 3 for g<i an d 53 observable. For each parameters, the averaged 
value is represented with a solid line whereas the dotted lines represent one standard deviations of 
the ensemble. It is clear that for K > 1 the parameters converge towards the theoretical values 
whereas for regular motions, the computed parameters are not representative of a GEV distribution 
and exhibit a spread that is more than five times bigger with respect to the chaotic counterpart. 
The results are similar if we change the initial conditions and the value of a. 

Let us now fix the value of K = 6.5. We want to show that we can depict the structure of the 
Standard map with the indicators presented above starting from 500x500 initial conditions uniformly 
distributed on the bidimensional torus. The number of iterations k, n, m and a are fixed as before. 
Results are shown in figure [2] where we compare the four parameters of GEV distribution (top and 
middle panels) with the Reversibility Error and Divergence of orbits in logarithm scale (lower panel). 
The number of iterations for the round off indicators is t = 100. It is evident that the structure 
of the Standard Map is well highlighted by all the indicators based on GEV distribution. For g\ 



the empirical values agreed with the theoretical ones £' = and a = 1/2 expected in the chaotic 
regions, whereas in the small regular islands we observe significant deviation from the expected 
values. Similar results hold for 52 and g% for which the expected theoretical values are £' = 1/6 
and £' = —1/6 respectively. The Round off based indicators highlight the same structure: in the 
chaotic sea the exponential growth of these quantities lead to a saturation at a value comparable 
with the size of the torus, whereas in the regular region their values remain order of magnitude lower. 

Eventually, we repeat the same analysis changing only the parameter K and fixing it to the value 
K = 4.5. The structure of the orbits in the standard map is now quite different, as we encounter 
larger stability islands and even the motions that start in the chaotic sea spend time around the 
regular islands. Due to this reason, we expect to observe a slower convergence to the parameter of 
the extreme value distribution with respect to the case presented above. This depends upon the 
fact that we may not extract a proper maximum if the motion is affected by a superposition of 
regular and chaotic motion. A better convergence must be observed if we increase the observations 
in each bin tom= 10 as an extreme value less affected by the regular motions can be extracted 
improving the convergence to the therotical parameters. This is exactly what is shown in figure [3] 
in the top panel £'(<7l) is compared for the case m = 10 3 (Left-hand side) and m = 10 4 (Right-hand 
side): clearly the convergence to the expected theoretical values £' = improves when the number 
of observations in each bin increases. Similar results hold for the others GEV parameters: as a 
further example £'(52) is presented in the middle panel of figure M The structure of the map is 
again well highlighted by GEV indicators if we compare them with the logarithm of Rt, At shown 
in the lower panel for t = 100. 

5 Conclusions 

Nowadays there exists a large family of dynamical indicators of stability that can be used to obtain 
informations about the stability of orbits in dynamical systems. Netherless, the numerical algorithm 
to compute many of them can be computationally expensive as it usually require the solution 
of variational equations or a Monte Carlo simulation. Moreover each indicator is specialized in 
detecting differences either in the chaotic regions or in the regular ones. 
In this framework we have introduced the parameters of the GEV distribution for some special class 



of observables to show that they effectively work as other dynamical indicators. As shown in Freitas 



et al. (2009), Faranda et al. (2011b) and Faranda et al. (2011a) a GEV like distribution can be 



observed only if the dynamical system satisfies certain mixing conditions or exhibit an exponential 
decay of the Hitting Time Statistics. 

We selected the parameters that do not depend explicitly on the number of observations in each bin 
m, so that our results only depend on the dimension of the attractor d. We proved the effectiveness 
of the indicators testing them on the Standard Map. First we proved that, varying the parameters 
K that regulate the chaoticity of the map, the parameters of the GEV fitted distribution for a small 
ensemble of initial conditions approach the theoretical values when K » 1 that corresponds to 
chaotic motions. We have also proved that GEV parameters are able to distinguish regular islands 
in chaotic sea for Standard Map parameter K = 6.5 and K = 4.5. In the latter case we have 
experienced the presence of superimposition of regular and chaotic motions that prevented us from 
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obtaining the theoretical expected parameters for the chaotic regions whem m was not big enough. 
To conclude, the GEV parameters provide basically the same informations of other indicators as the 
Divergence of orbits and the Reversibility error and they are easily accessible from a computational 
view point as they do not require the solution of the variational equations. With respect to the other 
indicators they highlight other relevant informations on the dynamics: the parameters are dependent 
on the dimension of the attractor that can be derived from a fit to the empirical distribution once 
we are sure that we are dealing with a chaotic orbit. Furthermore, the extreme value distribution 
is itself interesting from both practical and theoretical point of view: as far as the observable are 
concerned it gives detailed informations about the closest return near a certain initial condition. 
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Figure 1: Standard map: GEV parameters averaged over 500 different initial condition centered in 
x = 0.305, y = 0.7340 VS K. a) £'( 9l ), b) ?(g 2 ), c) £'( 53 ), d) a( 9l ). 
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Figure 3: Structure of the Standard Map for k = 4.5. a) Left-hand side: £'(<7i) for m,n = 10 3 , 
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